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Abstract 



Recently Burrage, de Rham, Heisenberg and ToUey have constructed eternal, classical 
solutions with closed timelike curves (CTCs) in a Galileon model coupled to an auxiliary 
scalar field. These theories contain at least two distinct metrics and, in configurations 
with CTCs, two distinct notions of locality. As usual, globally CTCs lead to pathologies 
including nonlocal constraints on the initial Cauchy data. Locally, with respect to the 
gravitational metric, we use a WKB approximation to explicitly construct small, short- 
wavelength perturbations without imposing the nonlocal constraints and observe that 
these perturbations do not grow and so do not lead to an instability. 
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1 The goal 



The central question in physics is: Given a configuration on an initial surface, local 
rules for evolving it to later times and boundary conditions at large spatial distances, 
what is the configuration at a future time? Of course in practice one can only know the 
initial conditions up to some bounded uncertainty. A theory is clearly only useful if this 
finite uncertainty in the initial conditions leads to a finite uncertainty in the results. In 
the case of most classical field theories, the rules for evolution are hyperbolic partial dif- 
ferential equations, and so it suffices that the initial surface be a spatial Cauchy surface. 
Generically small perturbations are described by the wave equation, and so small high 
wavenumber perturbations of the initial conditions lead to high frequency perturbations 
of the time evolution which, critically, have the same amplitude as the original pertur- 
bations and so remain finite. If on the other hand one tries to place initial data on a 
timelike surface and then evolve, small but arbitrarily high wavenumber perturbations in 
the initial conditions lead to arbitrarily large perturbations after a fixed amount of evolu- 
tion, therefore initial conditions placed on such surfaces cannot be used to determine the 
configuration elsewhere. 

In a universe with closed timelike curves (CTCs) there is an additional constraint. 
Evolving the surface into the future using local equations of motion, one can find the 
configurations on future surfaces. But iterating this process, one eventually returns to 
the initial surface itself. The physical observables must be single- valued, which means 
that the new values on the surface must be equal to the initial conditions. This constraint 
is in generally nonlocally encoded in the initial data, and so it presents a serious if not in- 
surmountable complication to the theory. One may be tempted to simply define an initial 
surface to which the CTCs are always parallel. Then these constraints can be satisfied by 
simply imposing a periodicity condition on the surface, such that circumnavigating the 
curve one arrives at the same values. However, as the curve is timelike, the surface is not 
spacelike and so again cannot be integrated to give the configuration elsewhere up to any 
finite error. If one tries to go ahead and try to integrate the equations of motion starting 
from initial data on such a surface, one will find a wrong sign dispersion relation leading 
to an apparently complex energy for the plane wave at high momenta, which indicates an 
instantly fatal instability. 

Does this apparent instability indicate that the configuration is really unstable, and so 
perhaps the CTC does not develop? Any configuration, even a free theory in Minkowski 
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space, exhibits such an instabihty if one attempts to integrate initial conditions defined 
on a timehke shce. Therefore, on its own, this feature does not imply an instability. 
On the contrary, if the equations are hyperbolic then a spacelike Cauchy surface can be 
constructed and its data integrated, leading to no obvious pathologies in the classical 
theory apart from the aforementioned nonlocal constraint. 

In this short note we would like to argue that the classical instability discovered in 
Ref. [I], in a variant of the Galileon model, is of this kind and so does not necessarily 
indicate that the classical theory cannot develop CTCs. While the authors did not directly 
impose initial conditions on a timelike slice, the periodicity condition which they use to 
solve the nonlocal constraints effectively fixes initial conditions for each mode. Physical 
high frequency instabilities occur when the timelike eigenvalue in the effective metric 
crosses zero and ghostlike instabilities when all eigenvalues change signs. On the other 
hand, in this case we will show that no eigenvalue of either relevant metric changes sign. 

The difference in points of view may appear to be only be semantic: On the one 
hand in Ref. [1], the authors observed that a small perturbation which obeys the nonlocal 
constraints leads to an instantaneous instability and so invalidates the solution. On the 
other hand in this note, using only spatial Cauchy surfaces, we will see in Sec. [3] that 
an infinitesimal deformation propagates with a constant amplitude, but most likely is 
inconsistent with the nonlocal constraints and so again invalidates the solution. However, 
there is a distinction between these two pathologies. In Ref. [T] it was claimed that the 
instability means that the theory avoids the CTC and so is consistent. The nonlocal 
constraint, on the other hand, yields a Galileon field which is generically is multiply 
defined and so the classical theory is inconsistent. Of course, if the Galileon is embedded 
in a healthy UV completion, this second problem is simply a sign that the low energy 
effective theory breaks down, but we will claim that this breakdown is not evident in 
any of the data within a local neighborhood where distances are determined using the 
nondynamical gravitational metric. Clearly an effective theory is more useful if there is a 
local criterion which determines when it can and cannot be trusted. 

2 History of CTCs in Galileon models 

The Galileon model is the most general theory of an interacting scalar field with a 
particular global symmetry under shifts and also linear transformations and with equa- 
tions of motion that only depend upon the first two derivatives of 0, guaranteeing the 
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absence of Ostrogradski ghosts. The Lagrangian density is 

1 

£ = ci£i - -df,(j)d''(f) + C3C3 + + cr^Cr, (2.1) 

where are real parameters and the Ck are given in Ref. [2]. For simphcity we will 
consider the special case Ci = C4 = C5 = which describes a limit of the 5-dimensional 
DGP model [3]. The only term Ck that we will then need is 

C3 = -\{nmA?- (2.2) 

The Galileon model is interesting in part because it contains ghost-free and classically 
stable solutions which violate the null energy condition [2] and even healthy solutions 
which cross from a domain in which they do not violate this condition into a domain 
in which they do |1] realizing the quintom scenario proposed for example in Ref. [S]. 
Recently it has been shown that if the Galileon symmetry is gauged, then the Galileon is 
the helicity zero component of a massive graviton [6], expanding upon the identification 
of the corresponding kinetic term in [7] . 

Both the DGP [8] and Galileon |9] models exhibit superluminal propagation. More 
precisely, small perturbations 

= 00 + (2.3) 

about a nontrivial background field configuration 0o are described by the wave equation 
in a background with an effective metric which differs from the metric of spacetime 

6C = -^(9^0)^(1 + 2c3n0o) + c3{d^d''<p^)d^6<pdM = -]^G^''d^5<pdM- (2.4) 

If 00 is harmonic with respect to the Minkowski gravitational metric then, normalizing 
so that C3 = 1, the inverse effective metric is just 

G^"" = 7]^"" - 9^9>o. (2.5) 

The scalar field propagates along null directions with respect to the effective metric, which 
in generic configurations can be timelike or spacelike with respect to the gravitational 
metric. Such superluminal propagation is not necessarily inconsistent so long as it does 
not lead to closed timelike curves [lOj, and in fact has recently become quite popular in 
phenomenology [H]. Nonetheless, it has been proposed in Ref. [12] that such superluminal 
propagation may be used to create configurations with CTCs, along the lines suggested 
in Ref. [13]. 

This proposal was followed in Ref. [H] where the authors constructed a classical 
Galileon solution, in the DGP case, which they claim is stable and develops CTCs. The 
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construction began from the observation that left-moving plane wave solutions 





h = f{x + t) (2.6) 

allow luminal propagation for right-moving Galileons and superluminal propagation for 
left-movers. More precisely, Galileons travelling in the same spatial direction as the plane 
wave travel along the direction proportional to (/" + l)t + (/" — l)x which is null with 
respect to the effective metric 

(2.7) 

but superluminal with respect to the gravitational metric if /" is negative. In fact, if 
/" < —1 then the Galileon travels backwards in time, which is a necessary condition at 
least somewhere in order to create a CTC. For simplicity the authors considered boundary 
conditions in which the Galileon and its first derivative vanish at both ends of straight 
rods, which means that /" cannot always be negative, as it must integrate to zero. The 
positivity of /" would lead to a massive subluminality which would eliminate the CTC. 
To avoid this problem, the authors included another plane wave moving to the right, 
for which the positive /" could be added without imposing subluminality on Galileons 
traveling to the left, since Galileons traveling in the opposite direction to the plane wave 
are always exactly luminal. The solution was thus complicated, but appeared consistent 
and stable. In fact, locally the effective metric was that of flat space pQ. 

In Ref. [T] the previous construction was simplifled by considering a single plane wave 
traveling around a compact circle. In this case, the well-deflnedness of the Galileon fleld 
clearly requires its second derivative to integrate to zero, as in the case of the previous 
paper. With a single plane wave, this leads to a very subluminal Galileon inside of the 
region in which the second derivative is positive. The authors avoided this problem by 
coupling the Galileon to an additional fleld x which is sensitive to the gravitational 
metric rj 



^G^''d,d(j)dM - \' 

where /i is a kinetic operator consisting of various contractions of derivatives, they argued 
that such couplings are generated in a variety of settings. The fleld x always travels 
luminally and so can bypass the region in which the Galileon fleld would be constrained 
to be subluminal without losing too much time to close the CTC. This CTC therefore 
can be circumnavigated only by two flelds which each feel different metrics, motivating 
their name "bimetric CTCs." 
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Generalizing the usual arguments of Hawking's chronology protection conjecture [15j 
the authors claimed that in the quantum theory, in a free field truncation, one expects a 
divergent backreaction which, depending on its sign, may inhibit the formation of CTCs. 
However they also claimed that the system would not arrive at this point, because already 
classically under arbitrarily small perturbations with wavenumber k it decays in a char- 
acteristic time 1/k, which for sufficiently high k is instantaneous. Of course if the CTCs 
do form, then this k is an initial condition on a temporal surface, and so the instability is 
just a result of the fact that the initial conditions were not defined on a spacelike surface. 
Before the formation of the curve there is no such instability but the energy measured 
with respect to the timelike isometry diverges. This should come as no surprise, as the 
corresponding Killing vector becomes lightlike as the CTCs (null Cauchy surface) form 
and so the corresponding frame is effectively infinitely boosted. In the next section we will 
review these solutions and determine their responses to such perturbations. We will see 
that there is always a choice of time direction such that these apparent local instabilities 
and divergences are absent from the classical theory. 

3 Bimetric CTCs 
3.1 Perturbations 

We will now review the bimetric CTCs of Ref. and study their classical stability 
with respect to small, local, high wavenumber perturbations. As we are interested in 
high wavenumber perturbations, a WKB approximation will be considered in which the 
magnitudes of the fields 6(f) and x change slowly, and so the differential operators act to 
leading order on their phases. More precisely, we will write these fields locally in the plane 
wave form 

with the approximation that the derivatives of the wavenumbers k are smaller than the 
various products fc^. 

Now the equations of motion are simply the algebraic condition that the vector ($, X) 
be annihilated by a matrix determined in terms of the functions k and /cj. Following pQ 
it is convenient to introduce the light cone momenta 

k± = k,± kt, p^ = kl + kl (3.2) 
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Then the algebraic equations of motion, to leading order in the WKB approximation, are 

+ k+k^ J \ X J ■ ^ ■ ^ 

The solutions to this equation correspond to values of the momenta for which the deter- 
minant of the matrix vanishes. This degeneracy condition can be solved for k^ 

k^-^ = ^ (-2p' + klf" ± ^Air + W) • (3.4) 

Recall that while /" on average is zero, it is less than —1 while the Galileon travels 
backwards in time, which is a necessary condition for the formation of a CTC Therefore 
there will be regions in which \f"\ > 1. The construction of Ref. [1] assumes that, at least 
inside of these regions 

4/i/i"^ « kl (3.5) 
and so in these regions one may approximate 

*- g it/"<o • - -1 i + k,r itr<o • '•^■^> 



Notice that, if /" > 0, then fci^"* solves (3.3) in the case X = fi = and k[_ ^ in the 
case $ = /i = 0. In other words, to leading order in /i^/A;^, the solution k^^^ describes a 
purely Galileon perturbation in the region in which /" > and a purely x perturbation 
in the region in which /" < 0. Recall that the Galileon is subluminal in the regime in 
which /" > 0, while x is always exactly luminal. Therefore the k^^^ fluctuations are never 
superluminal. This is the reason that no instability was found in Ref. fi\ when initial 
conditions were chosen for k^^^ on the CTC, because this curve is in fact spacelike with 
respect to the metric felt by k[^^ and so the Cauchy problem is well defined. On the other 
hand, k^S^ is nearly all Galileon in the regime /" < 0, where the Galileon is superluminal, 
and is nearly all x where the Galileon is subluminal. Therefore these are the modes which 
can circumnavigate the CTCs, and so those that were problematic in Ref. [1]. 

Between the regions in which /" is positive and negative, it of course must have a zero. 



In these regions the /i term dominates the square root in (3.4) and so Eq. (3.6) does not 
apply. In these regions, each solution interpolates between the Galileon dominated and 
X dominated regimes. This explains the observation in Ref. [1] that even for arbitrarily 
small /i, the behavior of these solutions does not reproduce that of the uncoupled case, 
in which no interpolation occurs. Any small, but non degenerate, kinetic operator will 
necessarily dominate in the square root near the zeros of /", leading to an interpolation 
which becomes significant as /" again becomes large. 
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3.2 Stability 



As fci^'* perturbations are not superluminal, and were shown to be healthy in Ref. [T], 
we will not consider them further. The authors found that, when fixing a wavenumber 
for k^~^ on the closed timelike curve, the evolution in another direction is unstable with 
respect to ky fluctuations. However, once the CTC has formed this curve is timelike for 
the k[_ ■* perturbations, and so fixing the wavenumber along the curve leaves an elliptic 
differential equation for the ky and perpendicular spatial directions. Elliptic differential 
equations cannot be integrated from a single initial slice due to just this instability. Such 
a pathology would be found if one imposed initial conditions on a timelike surface in any 
spacetime, and so on its own does not indicate a physical instability. In addition the 
divergent energy that they found as the closed timelike curves adiabatically form are seen 
in any background if one defines the energy with respect to a translation which is rotated 
from timelike to null. 

While a general search for instabilities is beyond the scope of this work, we will search 
for high frequency, classical, local instabilities. These depend only upon the local infor- 
mation available, in particular on the local value of the effective metric. High frequency 
instabilities occur when the effective metric becomes Euclidean, whereas ghost instabilities 
occur when it changes sign. We will show that no eigenvalues change sign. 



The two different mass-shell conditions in Eq. (3.6) indicate two different effective 
metrics, one in the Galileon dominated regime /" < and one in the x dominated 
regime /" > 0. By construction the x field is sensitive to the Minkowski metric, whose 
eigenvalues are constant and so no local instability is possible. Therefore we will focus on 
perturbations in the Galileon regime. 



The mass shell condition (3.6) indicates that the Galileon field is sensitive to the same 
inverse metric as in Ref. [H] 



fl + f" 
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-1 + f" 
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(3.7) 



The determinant of this metric is always equal to —1, for any value of /", therefore no 
eigenvalues can change sign and no instability occurs. 

The unstable modes discovered in Ref. [1] are waves traveling in the y direction. The 
authors imposed boundary conditions on the closed timelike curve, which quantized the k^ 
component. They found a high frequency instability for every value of k^-, so for simplicity 
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we will set = for the moment. In this case, a wave in the t — y plane feels the diagonal 
effective metric with values 1 + /" and —1. Clearly, when /" < —1, as it must be for the 
Galileon to travel backwards in time, this metric becomes Euclidean and a high frequency 
instability occurs. Equivalently, the energy computed using the derivative with respect 
to the t direction is complex when +/" < — 1 and diverges as +/" approaches —1. 

However, this instability relies critically on fixing k^, which corresponds to fixing an 
initial condition on the timelike surface. In fact, the x direction is timelike precisely when 
/" < —1. What happens if one does not fix k^7 Is there a stable oscillation mode for 
any fixed value of kyl In other words, if one perturbs the Galileon field locally with a 
wave which rapidly oscillates in the y direction and lets it evolve, what happens to the 
amplitude of the oscillation? Of course, to define this question, one must define evolve. 
Locally, evolve means to integrate the equations of motion in a temporal direction, that 
is a direction whose norm with respect to the effective metric is positive. In our WKB 
approximation this is equivalent to solving the mass shell conditions for u, where u is the 
frequency with respect to a positive-normed direction. 

One direction whose norm is certainly positive is the timelike eigenvector of the inverse 
effective metric 



{h,k,) = {f",-l + Vf"" + l)- (3.8) 
Therefore the mass-shell condition 

= k^G^'k, = (1 + f")k^ + 2f"ktk, + (-1 + f")kl + kl + kl (3.9) 

may be satisfied for a high frequency perturbation which oscillates in the y direction with 
an arbitrary momentum ky by simply adding to k a. multiple of this timelike eigenvector 



-Kf" ky{i-^r^i) 

:, Abj^, U 



2(/"2 - /" + l)(/" + /pTl) ^2(/"2 _ f" + l)(/" + /pTI) 

(3.10) 

This wavevector, in the WKB approximation (3.1), describes the classical evolution of 
a local, high-frequency Galileon perturbation which oscillates in the y direction with 
arbitrary wavenumber ky. The positivity of the expressions in the various square roots of 



(3.10) imply that all components of k are real, and so no local instability is present for 
any value of /", at least to leading order in the WKB approximation. 

Of course the vanishing proper length of this curve at the moment at which it becomes 
null means that another notion of locality exists for the Galileon field, using the effective 
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metric, in which the entire curve is locaQ Instabihties which are local in that sense are 
considerably more difficult to analyze, resembling the situation with CTCs in gravity. 
As described in Refs. [151 lini H] the vanishing proper distance may indicate a divergent 
stress tensor in the quantum theory, which certainly would lead to a large backreaction 
in the gravitational theory and perhaps also in this nongravitational theory. However the 
sign of this backreaction is important. It may be that it prohibits the formation of the 
CTC as was speculated in the case of free fields in Refs [13, ttZl [T] , or it may be that the 
UV theory intervenes only to cut off the divergence without inhibiting the formation of 
CTCs [IE]. 
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